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Abstract. We study the computational question whether a given polytope or spec- 
trahedron Sa (as given by the positive semidefiniteness region of a linear matrix pencil 
A(x)) is contained in another one Sb- 

First we classify the computational complexity, extending results on the polytope/poly- 
tope-case by Gritzmann and Klee to the polytope/spectrahedron-case. For various re- 
stricted containment problems, NP-hardness is shown. 

We then study in detail semidefinite conditions to certify containment, building upon 
work by Bcn-Tal, Nemirovski and Helton, Klep, McCullough. In particular, we discuss 
variations of a sufficient semidefinite condition to certify containment of a spectrahe- 
dron in a spectrahedron. It is shown that these sufficient conditions even provide exact 
semidefinite characterizations for containment in several important cases, including con- 
tainment of a spectrahedron in a polyhedron. Moreover, in the case of bounded Sa the 
criteria will always succeed in certifying containment of some scaled spectrahedron vSa 
in S B . 



1. Introduction 

Denote by Sk the set of all real symmetric k x /c-matrices and by Sk[x] the set of 
symmetric k x k- matrices with polynomial entries in x = {x\, . . . , x n ). For Aq, . . . , A n G Sk, 
let A{x) denote the linear (matrix) pencil A(x) = Aq + x\Ai + • • - + x n A n G <Sfe[x]. Then 
the set 

(1.1) S A := {x G R n : A{x) h 0} 

is called a spectrahedron, where A(x) >z denotes positive semidefiniteness of the matrix 
A(x). 

Spectrahedra arise as feasible sets of semidefinite programming (see [HI ES] ) - In the last 
years, there has been strong interest in understanding the geometry of spectrahedra (see, 
e.g., [H QUI EE] ) , particularly driven by their intrinsic relevance in polynomial optimization 
[ilfTT] and convex algebraic geometry [19], [20] . Spectrahedra naturally generalize the class 
of polyhedra, see [3j [27] for particular connections between these two classes. 

In this paper, we study containment problems for polyhedra and spectrahedra. Since 
polyhedra are special cases of spectrahedra, we can use the following general setup: Given 
two linear pencils A(x) G Sk[x] and B{x) G Si[x], is Sa Q Sb^ 

For polytopes (i.e., bounded polyhedra), the computational geometry and computa- 
tional complexity of containment problems have been studied in detail. See in particular 
the classifications by Gritzmann and Klee fT2 \ [T31 [Tlj . Notably, it is well-known that the 
computational complexity of deciding containment problems strongly depends on the type 
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of the input. For instance, if both polytopes are given by their vertices (V-polytopes), or 
both polytopes are given as an intersection of halfspaces ("H-polytopes), containment can 
be decided in polynomial time, while it is co-NP-hard to decide whether an "H-polytope 
is contained in a V-polytope (see [HE])- 

For spectrahedra, much less is known. Ben-Tal and Nemirovski studied the matrix cube 
problem [2J, which corresponds to the containment problem where Sa is a cube. In a much 
more general setting, Helton, Klep, and McCullough [16J studied containment problems 
of matricial positivity domains (which live in a union of spaces of different dimensions). 
As a byproduct, they also derive some implications for containment of spectrahedra. 

In the current paper, we study containment problems of polytopes, polyhedra and 
spectrahedra from a computational viewpoint. In Section [31 we extend existing com- 
plexity classifications for the polyhedral situation to the situation where polytopes and 
spectrahedra are involved. In particular, the containment question of a V-polytope in a 
spectrahedron can be decided in polynomial time, and the question whether a spectrahe- 
dron is contained in an "H-polytope can be formulated by the complement of semidefinite 
feasibility problems (involving also strict inequalities). Roughly speaking, the other cases 
are co-NP-hard. This includes the containment problem of an H-polytope in a spectra- 
hedron, already when the spectrahedron is a ball. The complete classification is stated in 
Theorems EOlfel 

To overcome the situation that the general containment problem for spectrahedra is 
co-NP-hard, relaxation techniques are of particular interest. Our point of departure in 
Section H] is the relaxation from [16] which provides a distinguished sufficient criterion 
for containment of a spectrahedron Sa in a spectrahedron Sb (see (14. 2p ). We provide an 
elementary derivation of this semidefinite relaxation (as opposed to the operator-theoretic 
techniques used there) and study the quality of the criterion. This leads to a new and 
systematic access to studying containment problems of polyhedra and spectrahedra and 
provides several new (and partially unexpected) results. 

In particular, we obtain the following new results: 

1. We discuss variations of containment criterion (14. 2J , which lead to improved numer- 
ical results, see Theorem 14.31 Corollaries 14.61 and 14.71 and Section 16 . 11 

2. We exhibit several cases when the criteria are exact (see Theorem 14. 8p . For some of 
the cases we can provide elementary proofs. The main case in Theorem 14.81 states that the 
sufficient criteria for the containment of spectrahedra in polyhedra (in normal form) are 
exact characterizations. The proof of the statements is given in Section by developing 
various properties of the containment criteria (transitivity, block diagonalization) and 
combining them with duality theory of semidefinite programming. The exactness of the 
spectrahedron-polyhedron-case is particularly surprising, since a priori the criteria depend 
on the linear pencil representation of the spectrahedron. 

3. In Section I6.1[ we extend the results from [16] on cases, where the criteria are not 
exact. For a counterexample in [TB] we exhibit the phenomenon that the containment 
criteria will at least succeed in certifying that a scaled version of the spectrahedron Sa is 
contained in Sb- 
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4. In Proposition I6.2[ we show that in the case of bounded Sa there always exists a 
scaling factor v > such that for the scaled spectrahedron pair (uSa, Sb) the criteria (14.21) 
and g3D hold. 

We will close the paper by explaining some implications of the scaling result on the 
optimization version of containment problems (as also relevant for the computation of 
geometric radii of convex bodies, e.g., in [T21 US])- 



2. Preliminaries 

Throughout the paper we work in n-dimensional Euclidean space M n , and || • || denotes 
the Euclidean norm. 



Matrices and block matrices. For a matrix A, the (i, j)-th entry of A is labeled by a^ 
as usual. For a block matrix B, we label the (i,j)-th block by and the (s, £)-th entry 
of Bij by (B^st. 

A square matrix with 1 in the entry and zeros otherwise is denoted by E%j. The 
n x n identity matrix is denoted by I n . 

The Kronecker product A <g> B of square matrices A of size k x k and B of size I x I is 
the kl x kl matrix 



(2-1) 



A® B 



'an B 



Oifc B' 



dkk B_ 



(see, e.g., [H1|2T]). It is well-known (see, e.g., [2TJ Cor. 4.2.13]) that the Kronecker product 
of two positive semidefmite matrices is again positive semidefinite. 



Polyhedra and polytopes. A polyhedron is the intersection of finitely many halfspaces. 
A bounded polyhedron or, equivalently, the convex hull of finitely many points in R n is 
called polytope. 

For algorithmic questions in n-dimensional space it is crucial whether a polytope is 
given in the first way {H-polytope) or in the second way (V-polytope). Our model of 
computation is the binary Turing machine: polytopes are presented by certain rational 
numbers, and the size of the input is defined as the length of the binary encoding of the 
input data (see, e.g., [12]). A V-polytope P is given by a tuple (n; m; . . . , i/ m )) with 
n, m G N, and . . . , v^ m > G Q n such that P = conv{v^\ . . . , v^ m '}. An "H-polytope P 
is given by a tuple (n; m; A; b) with n, m G N, a rational m x n-matrix A, and b G Q m 
such that P = {x G M. n : b + Ax > 0} is bounded. If the i-th row (b + Ax)i > defines a 
facet of P, then the z-th row of A is an inner normal vector of this facet. 

For fixed dimension, "H- and V-presentations of a rational polytope can be converted 
into each other in polynomial time. In general dimension (i.e., if the dimension is not 
fixed but part of the input) the size of one presentation can be exponential in the size of 
the other [24]. 
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Spectrahedra. Given a linear pencil 

■n 

(2.2) A(x) = A + Y^ x p A p e S k [x) with A p = (a?-) , < p < n , 

P =i 

the spectrahedron = {i 6 R n : A(x) ^0} contains the origin in its interior if and only 
if there is another linear pencil A'(x) with the same positivity domain such that A' = I, 
see [51 ED]- In particular, then Sa is full-dimensional. To simplify notation, we sometimes 
assume that A(x) is a monic linear pencil, i.e. Aq = I k . As a shorthand we use A y B to 
state that A — B is positive semidefmite. 

Note that every polyhedron P = {i 6 1™ : b + Ax > 0} has a natural representation 
as a spectrahedron: 

( rai(x) " 

(2.3) P = P A = lxeR n : = ' • • ^0 

I LoO o fe (x). 

where aj(x) abbreviates the i-th entry of the vector 6 + Ac. Pa contains the origin if 
and only if the inequalities can be scaled so that b = t k , where 1^ denotes the all-ones 
vector in M fc . Hence, in this case, A(x) is monic, and it is called the normal form of the 
polyhedron Pa- 

A centrally-symmetric ellipsoid with axis- aligned semi- axes of lengths a%, . . . , a n can be 
written as the spectrahedron Sa of the monic linear pencil 



(2.4) A(x) = I n+1 + J2—( E P,n+i + E n+l , p ). 

ttp 
p=i y 

We call (12 .4p the normal form of the ellipsoid. Specifically, for the case of equal semi-axis 
lengths r := ai = ■ ■ ■ = a n this gives the normal form of a ball with radius r. 

For algorithmic questions, a linear pencil is given by a tuple (n; k; Aq, . . . , A n ) with 
n, k e N and Aq, . . . , A n rational symmetric matrices. 

3. Complexity of containment problems for spectrahedra 

In this section, we classify the complexity of several natural containment problems 
for spectrahedra. For polytopes the computational complexity of containment problems 
strongly depends on the type of input representations. For V- and ^-presented polytopes, 
the following result is well-known (see [7f [14"]). 

Proposition 3.1. Deciding whether a polytope P is contained in a polytope Q can be done 
in polynomial time for the following cases: 

(1) Both P and Q are 'K-polytopes , 

(2) both P and Q are V-polytopes, or 

(3) P is a V -polytope while Q is an H-polytope. 

However, deciding whether an %-polytope is contained in a V -polytope is co-NP- complete. 
This hardness persists if P is restricted to be a standard cube and Q is restricted to be the 
affine image of a cross polytope. 
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In the next statements we extend this classification to containment problems involving 
polytopes and spectrahedra. See Table [I] for a summary. Theorems 13.21 and 13.31 give the 
positive results. 







V 


S 


n 


p 


co-NP-complete 


co-NP-hard 


V 


p 


P 


P 


s 


"SDP" 


co-NP-hard 


co-NP-hard 



Table 1. Computational complexity of containment problems, where the 
rows refer to the inner set and the columns to the outer set and S abbreviates 
spectrahedron. 



Theorem 3.2. Deciding whether a V-polytope is contained in a spectrahedron can be done 
in polynomial time. 

Proof. Given a V-presentation P = conv{?;W, . . . ,v^ m '} and a linear matrix pencil A(x), 
we have P C Sa if and only if all the points are contained in Sa- Thus, the containment 
problem is reduced to m tests whether a certain rational matrix is positive semidefinite. 
This can be decided in polynomial time, as one can compute, for a rational, symmetric 
matrix A, a decomposition A = UDU T with a diagonal matrix D in polynomial time (see, 
e.g., 0). ' □ 

Containment questions for spectrahedra are connected to feasibility questions of semi- 
definite programs in a natural way. A Semidefinite Feasibility Problem (SDFP) is defined 
as the following decision problem (see, e.g., [26]): Given a linear pencil defined by a tuple 
(n; k; Aq, . . . , A n ) with n, k G N and Aq, . . . ,A n rational symmetric matrices. Are there 
real numbers xi,...,x n such that A(x) = A + Ylp=i x pA P h 0, or equivalently, is the 
spectrahedron Sa non-empty? 

Although semidefinite programs can be approximated up to an additive error of e 
in polynomial time, the question "SDFP G P?" is one of the major open complexity 
questions in semidefinite programming (see [6j,[26]). Consequently, the following statement 
on containment of a spectrahedron in an "H-polytope does not give a complete answer 
concerning polynomial solvability of these containment questions in the Turing machine 
model. If the additional inequalities were non-strict, then we had to decide a finite set of 
problems from the complement of the class SDFP. 

Theorem 3.3. The problem of deciding whether a spectrahedron is contained in an %- 
polytope can be formulated by the complement of semidefinite feasibility problems (involving 
also strict inequalities), whose sizes are polynomial in the description size of the input data. 

Proof. Let A(x) be a linear matrix pencil and P = {x G M n : b + Bx > 0} with 
B G Q mxn be an "H-polytope. For each i G {1, . . . , m} incorporate the linear condition 
bi + YTj=i^ij x j < i 11 ^ the linear pencil A(x). If one of the resulting m ("semi-open") 
spectrahedra is nonempty then Sa 2 P- D 
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The positive results in Theorems 13.21 and 13.31 are contrasted by the following hardness 
results. 

Theorem 3.4. 

(1) Deciding whether a spectrahedron is contained in a V-polytope is co-NP-hard. 

(2) Deciding whether an H-polytope or a spectrahedron is contained in a spectrahedron 
is co-NP-hard. This hardness statement persists if the H-polytope is a standard 
cube or if the outer spectrahedron is a ball. 

Proof. Deciding whether a spectrahedron Sa is contained in a V-polytope is co-NP-hard 
since already deciding whether an %-polytope is contained in a V-polytope is co-NP-hard 
by Proposition 13.11 

Concerning the second statement, co-NP-hardness of containment of H-polytopes in 
spectrahedra follows from Ben-Tal and Nemirovski [2J Proposition 4.1], who use a re- 
duction from the maximization of a positive semidefmite quadratic form over the unit 
cube. 

For the co-NP-hardness of containment of an "H-polytope in a ball, we provide a re- 
duction from the NP-complete 3-satisfiability problem (3-SAT |5J): Does a given Boolean 
formula <3> over the variables z\, . . . , z n in conjunctive normal form, where each clause has 
at most 3 literals, admit an assignment that evaluates True? 

The 2 n possible assignments {False, True}" for be identified with the 

vertices of an n- dimensional cube [—1,1]™. Let B be a ball (which is a spectrahedron), 
such that the vertices of [—1, l] n just "peak" through its boundary sphere S. Precisely 
(assuming w.l.o.g. n > 2), choose the radius r of B such that 



Note that such a radius can be determined in polynomial time and size. 

For the definition of the "H-polytope P, we start from the %-presentation {x G M. n : 
— 1 < Xi < 1, 1 < i < n} of [—1, l] n and add one inequality for each clause of Let 
C = C\ V ■ ■ ■ VC m be a 3-SAT formula with clauses C\ , . . . , C m . Denote by zl the complement 
of a variable Zj, and define the literals z\ := Zi, z® := ~zi. If the clause Ci is of the form 
Ci = z^ 1 V z^ 2 V z^ with ej i; ej 2 , G {0, 1} then add the inequality 



If P C B then, by the choice of r, none of the points in {—1, 1}™ can be contained in 
P and thus there does not exist a valid assignment for $. Conversely, assume that P is 
not contained in B. Let p G P \ B C [— 1, l] n . We claim that componentwise rounding 
of p yields an integer point p' G { — 1, l} n satisfying all defining inequalities of P. To see 
this, first note that by the choice of the radius of B, the components pi of p differ at most 
e < < | from either —1 or 1. 

In order to inspect what happens to the inequalities when rounding, assume without 
loss of generality that the inequality is of the form x\ + x 2 + x 3 > —1. We assume a 
rounded vector p' does not satisfy the inequality, even though p does: 




(-l) e ^X h +{-l)^X i2 + (-l) ei sx i3 < 1. 



(3.1) 



Pi +P2+P3 < - 1 > but P1+P2+P3 > - 1 • 
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Since p' G {—1, l} n , ( 13. ip implies p[ = p' 2 = p' 3 = — 1. Hence, at least one of p±,p2 and 
p% differs from either —1 or 1 by more than 1/3, which is a contradiction. This completes 
the reduction from 3-SAT. 

Finally, deciding whether a spectrahedron Sa is contained in a spectrahedron Sb is co- 
NP-hard, since already deciding whether an "H-polytope is contained in a spectrahedron 
is co-NP-hard. □ 

4. Relaxations and exact cases 

In this section, we revisit and extend the relaxation techniques for the containment of 
spectrahedra from [16]. Our point of departure is the containment problem for pairs of 
■H-polytopes, which by Proposition 13.11 can be decided in polynomial time. Indeed, this 
can be achieved by solving a linear program, as reviewed by the following necessary and 
sufficient criterion. 

Proposition 4.1. Let P A = {x G R n : t k + Ax > 0} and P B = {x G R n : 1; + Bx > 0} 

be polytopes. There exists a right stochastic matrix C (nonnegative entries, each row 
summing to one) with B = CA if and only if Pa Q Pb- 

For preparing related statements in more general contexts below, we review the proof 
which uses the following version of Farkas' Lemma: 

Proposition 4.2 (Affine form of Farkas' Lemma [28j Corollary 7.1h]). Let the polyhedron 
P = {x G R™ : k(x) > 0, % = 1, . . . , m} with affine functions U : M n — > R be nonempty. 
Then every affine I : R n — > R that is nonnegative on P can be written as l(x) = cq + 
YliLi c ih{%) with nonnegative coefficients C{. 

Proof, (of Proposition 14.11 ) If B = CA with a right stochastic matrix C, then for any 
x G Pa we have 1; + Bx = t t + C(Ax) > 0, i.e., P A C P B . 

Conversely, if Pa Q Pb, then for any i G {1, . . . , 1} the i-th row (1; + Bx\ of 1; + Bx 
is nonnegative on Pa- Hence, by Proposition 14.21 (1; + Bx\ can be written as a linear 
combination 

k 

(li + Bx)i = l + {Bx)i = c' iQ + J2 c 'i j ^k + Ax) J 

3=1 

with nonnegative coefficients c'^. Comparing coefficients yields Y^!j=i c 'ij = c 'io- Since 
Pa is a polytope with zero in its interior, the vertices of the polar polytope P% are given 
by the rows — Aj of —A. Hence, for every i G {1, ...,/} there exists a convex combination 
= Xlj=i Aj) with nonnegative Ajj and Y^=t ^ij = 1) which we write as an identity 
SjLi^ij(lfc + Ax)j = 1 of affine functions. By multiplying that equation with c' i0 , we 
obtain nonnegative c'-j with Ylj=i c ij(^-k + Ax)j = c' i0 , which yields 

k 

l + (Bx)i = ^(^ + cg.)(l Jk + Ar) i . 
i=i 

Hence, C = (qj) with c^- := c\ • + cf - is a right stochastic matrix with B = CA. □ 
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The sufficiency part from Proposition 14. II can be extended to the case of spectrahedra in 
a natural way. The natural description of a polytope P as a spectrahedron, as introduced 
in Section [21 is given by 

p = p A = { x Gl" : A(x) = diag(ai(a;), . . . , a k {x)) y 0} , 

where a»(x) is the i-th entry of the vector + Ax. Then, as in the definition of a linear 
pencil (12. 2p . A p is the k x k diagonal matrix diag(A jP ) of the p-th column of A. Proceed 
in the same way with Pg. Now define a kl x kl matrix C by writing the entries of C on 
the diagonal, i.e. C = diag(cn, . . . , Qi, Ci 2 , . . . , q 2 , . . . , c^, . . . , c^). Then the condition 
from Proposition 14.11 translates to 

k k 

(4.1) C diagonal, C = (<%) J. =1 b 0, Ji = £ = 1, . . . , n : B p = £ , 

i=l i=l 

where af- is the (i,j)-th entry of A p and C?. G R' xZ . Theorem 14.31 below tells us, that 
C" does not need to be diagonal and yields a sufficient condition for the containment of 
spectrahedra. 

4.1. A sufficient condition for containment of a spectrahedron in a spectrahe- 
dron. Let A(x) G S]c[x] and B{x) G Si[x] be linear pencils. 

In the following, the indeterminate matrix C = {Cij) i - =l ("Choi matrix") is a sym- 
metric kl x ^/-matrix where the Cy are I x /-blocks. By showing the equivalence of 
containment of the so-called matricial relaxations of two spectrahedra Sa, Sb given 
by monic linear pencils and the existence of a completely positive unital linear map 
r : span{A , Ai, . . . , A n } — > span{£> , B>i, . . . , B n }, A p \-t B p , the authors of [151 US] 
proved that the system 

k 

(4.2) c = (C,,)': , , yo, \/p = 0, . . . , n ■ B p = J2 <A 

has a solution if and only if the matricial relaxation of Sa is contained in the one of Sb- If 
so, then Sa Q Sb- We show the latter statement in an elementary way, see Theorem 14.31 
Moreover, in our approach it becomes apparent that we can relax the criterion given by 
Helton, Klep and McCullough by replacing the linear constraint on the constant matrices 
in (14. 2 p with semidefinite constraints, 

k k 

(4.3) C=(C4 =1 b0, Bb-XXCktO, V /' 1 " : 'h- 

i,j=l i,j=l 

This relaxed system is still sufficient for containment of spectrahedra as the following 
theorem shows. 

Theorem 4.3. Let A(x) G Sk[x] and B(x) G Si[x] be linear pencils. If one of the sys- 
tems (14. 2p or (14. 3 p is feasible then Sa Q Sb- 
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Proof. First we show the statement for (14.31) . 

For x G Sa, the last two conditions in (14. 3 P imply 

n k n k k 

(4.4) B( X ) = b + e xp b p h E a l C * + E E ^ °5 ^ = E w*))** c « • 

p=l i,j=l p=l i,j=\ i,j=l 

Since A(x) and C are positive semidefinite, the Kronecker product A(x) ® C is positive 
semidefinite as well, see (12. ip . As a consequence, all principal submatrices of A(x) <8> C 
are positive semidefinite. Consider the principal submatrix where we take the (z,j)-th 
sub-block of every (i, j)-th block, 

((A(a:))yCy) fc G5„N. 

V Aj=l 

To be more precise, A(x) <8> C is a A; 2 / x A; 2 /-matrix with k x k blocks of the form 

~(A(x))ijC n ■ ■ ■ {A(x)) ij C lk ' 

(A(x)) ij C= ! (.l(.r)),/', ; ! e5 w . 

_(A(x)) y C fcl • ■ ■ (A(x))ijC kK 

(Remember that (A(x))ij is a scalar). Now we take the (i,j)-th block of (A(x))ijC, i.e. 
[A[x))ijCij. 

Set I = [Ii, . . . , Ii] T . Then the claim for system (14. 3p follows from the fact that positive 
semidefiniteness 'V " is a transitive relation on the space of symmetric matrices, that is, 

(4.5) v T B(x)v > v T (l T ((A(a;)) y C ^.._ x A v 



(v T ... ^)((^)) y .^) fc (« ... t;) T >0 

V / 1,7=1 



for every !)6l'. 

Specializing 'V" to "=" in (14. 4 p and ">" to "=" in (14. 5 p provides a streamlined proof 
for (Ojl . □ 

For both systems (14. 2 p and (14. 3 p the feasibility depends on the linear pencil represen- 
tation of the sets involved. In Section 16711 we will take a closer look at this fact. 

Remark 4.4. The sub-block argument in the proof of Theorem 14.31 can also be stated 
in terms of the Khatri-Rao product (see [23]). Let A = (Aij)ij and B = (-Bjj)jj be block 
matrices, consisting of k x k blocks of size p x p and q x q, respectively. The Khatri-Rao 
product of A and B is defined as the blockwise Kronecker product of A and B, i.e., 

A* B = { : \ jj ® Ih i ), i • 

If both A and B are positive semidefinite, then the Khatri-Rao product A * B is positive 
semidefinite as well, see [23, Theorem 5]. 

Now consider A(x) and C as in the proof of Theorem 14.31 Then p — 1 and q = I. 
Therefore, 

A(x) * C = ((A(x)) l3 ® C^l , = ((^(z))yQj)*j =1 
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is positive semidefinite. 

The subsequent statement shows that Theorem 14.31 is invariant under translation. Let 
Sa be a spectrahedron defined by the linear pencil A(x) = A + J2 p =i x pA p . To move Sa 
by a vector v = (f i, . . . , v n ) we substitute x — v into the pencil and get 

n n 

A(x -v) = A -y j v p A p + XpAp. 

p=i p=i 

Lemma 4.5 (Translation symmetry). The criteria (14. 2p and (14. 3p are invariant under 
translation. 

Proof. Given linear pencils A(x) and B(x), let C be a solution to system ( 14. 3p . Then it 
is also a solution for the translated pencils A(x — v) and B{x — v ) for any v £ M. n . The 
translation only has an impact on the constant matrix, we have to show 

n / k / n \ \ 

(4.6) b -J2 v p b p - (J2 K - E v$ °« ^ °- 

p=i \tj=i V p=i / / 

Since 5 P = ^ij=i a fjC*i f° r all p = 1, . . . , n, (14.61) is equivalent to B - J2i,j=i a %^ij ^ 0> 
which is the condition on the constant matrices before translating. 

As in the proof of Theorem 14.31 specializing u y" to "=" yields a proof for (14. 2p . □ 

If Sb is contained in the positive orthant, we can give a stronger version of the criterion 
introduced in Theorem 14.31 

Corollary 4.6. Let A(x) £ Sk[x] and B(x) £ Si[x] be linear pencils and let Sa be con- 
tained in the positive orthant. If the following system is feasible then Sa ^ Sb- 

k k 

(4.7) C = (C^ J=1 hO, B - J2 <A b 0, Vp = 1, . . . , n : B p - <fa t 0. 

i,j=l i,j=l 

Proof. The proof is along the lines of the proof of Theorem 14.31 Indeed, since Sa lies in 
the positive orthant, we have x > for all x £ Sa and hence, 

n k n k k 

B(x) = B + J2 x P Bp > £ a% Cij + £ £ x p a% C y = £ C y . 

p=l tj=l p=l i,j=l i,j=l 

□ 

By relaxing system ( 14. 2 p to ( 14. 7p the number of scalar variables remains ~kl(kl + 1), 
whereas the |(n+l)/(/ + l) linear constraints are replaced by n + 1 semidefinite constraints 
of size I x I. 

If containment restricted to the positive orthant implies containment everywhere else, 
criterion (14. 7ft can be applied, even if the spectrahedron is not completely contained in 
the positive orthant. To make use of this fact, we have to premise a certain structure of 
the spectrahedra. We give an example in the following corollary. 
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Corollary 4.7. Let A(x) G Sk[x] and B(x) G Si[x] be linear pencils defining spectrahedra 
with a reflection symmetry with respect to all coordinate hyperplanes. If system ( 14. 7p is 
feasible then Sa Q Sb- 

In Section 16.11 we will see that the relaxed version ( 14. 7p is strictly stronger than sys- 
tem ( 14. 2p . There are cases, where a solution to the relaxed problem ( 14. 7p exists, even 
though the original problem ( 14. 2 p is infeasible. 

4.2. Exact cases. It turns out that the sufficient semidefmite criteria ( 14. 2 p and (14. 3 p 
even provide exact containment characterizations in several important cases. Detailed 
statements of these results and their proofs will be given in Statements I4.9[ I4.11[ I4.13[ 

EU 

For ease of notation, most statements in this section as well as in section are given for 
monic pencils and proved only for criterion ( 14. 2p . Note however, that feasibility of (14. 2p 
implies feasibility of (14. 3p . Furthermore, after translating the (in this section mostly 
bounded) spectrahedra to the positive orthant, Corollary 14.61 can be applied. Since crite- 
rion ( 14. 2 p is invariant under translation, its feasibility again implies that system (14 .7p has 
a solution for the translated spectrahedra. 

Besides the normal forms of polyhedra, ellipsoids, and balls introduced in Section [2J the 
exact characterizations will also use the following extended form of a spectrahedron Sa- 
Given a linear pencil A(x) G «Sfc[x], we call the linear pencil with an additional 1 on the 
diagonal 

1 

A{x) 

the extended linear pencil of Sa = S^ (the spectrahedra coincide, only the representations 
of Sa and S^ differ, since the 1 we add for technical reasons is redundant). The entries of 
A p in the pencil A(x) = A + Ylp=i x pA p are denoted by af- for i,j = 0,...,k,as usual. 

Theorem 4.8. Let A(x) G Sk[x] and B(x) G Si[x] be monic linear pencils. In the 
following cases the criteria ( 14. 2 p as well as ( 14. 3 p are necessary and sufficient for the 
inclusion Sa Q Sb-' 

(1) if A(x) and B(x) are normal forms of ellipsoids (both centrally symmetric, axis- 
aligned semi-axes), 

(2) if A(x) and B(x) are normal forms of a ball and an %- polyhedron, respectively, 

(3) if B(x) is the normal form of a poly tope, 

(4) if A(x) is the extended form of a spectrahedron and B(x) is the normal form of a 
polyhedron. 

In this section, we provide the proofs of (1), (2), where the sufficiency parts follow 
by Theorem 14.31 The cases (3) and (4) will be treated in Section [5j We start with the 
containment of "H-polyhedra in "H-polyhedra which slightly generalizes Proposition 14.11 
and will be used in the proofs of later statements. 

Lemma 4.9. Let A(x) G Sk[x] be the normal form of a polyhedron as defined in (12. 3p 
and let A(x) = diag(l, ai(x), . . . , ak(x)) G be the extended linear pencil of A(x). 



(4.8) A{x) 



G S k+ i[x] 
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Let B(x) G Si[x] be the normal form of a polyhedron. When applied to the pencils A(x) 
and B(x) criterion (14. 2 j) is necessary and sufficient for the inclusion = Sa Q Sb- If 
Sa is a polytope, i.e. a bounded polyhedron, the pencil A(x) can be used instead of A(x). 

Proof. With regard to f 12 .31) . the polyhedra Sa and Sb are of the form Sa = {% G W 1 : 
t k + Ax > 0} and S B = {x G R n : h + Bx > 0}, respectively, and let A be the (k + l)xn 
matrix defined by A : = ^ 

If S% C 5s, then there exist convex combinations (ti+Bx)i = Qo+X]j=i c ij(^-k+Ax)j = 

Y^=o c *i(lfc+i + where (Ijt+i + Ax)o = 1 and the coefficients Cjj are nonnegative, 

just as in the proof of Proposition 14.11 

Now we construct a matrix C that is a solution to system ( 14.21) . Recall that C consists 
of matrices of size I x I : C = (C s t)g t=Q . Set the z-th diagonal entry of Cjj to be c^, 
and choose all other entries to be zero. The resulting matrix is a diagonal matrix with 
non-negative entries, which makes it positive semidefinite. 

Clearly, B(x) = Yl^o^k+i + A x)jCjj = Z)ij=o(lfc+i + Ax)jC { j. Comparing coeffi- 
cients, we see that Ii = Ylj=o Cjj and B p = J2 i j =0 afjCij for all p G {1, . . . , n}. 

If the inner polyhedron Sa is a polytope, the constant 1 is a convex combination of the 
remaining polynomials ai(x), . . . , ak(x). Thus the additional 1 in the upper left entry of 
pencil A(x) is not needed. □ 

As we have seen in the proof of Lemma H~9l for polyhedra there is a diagonal solution to 
(14.21) . Thus it is sufficient to check the feasibility of the restriction of (14.21) to the diagonal 
and checking inclusion of polyhedra reduces to a linear program. 

Remark 4.10. For unbounded polyhedra, the extended normal form is required in order 
for the criterion to be exact. Without it, already in the simple case of two half spaces 
defined by two parallel hyperplanes, system (14. 2 j) is not feasible. 

The following statement on ellipsoids uses the normal form (12. 4p . 

Lemma 4.11. Let two ellipsoids Sa and Sb be centered at the origin with semi-axes 
parallel to the coordinate axes, given by the normal forms 

n n 

A(x) = I n+1 + ^2 —{E P ,n+i + E n+1>p ) and B(x) = I n+1 + -^-(E Pin+1 + E n+1)P ) , 
P =i a P P =i b P 

respectively. Here (ai, . . . , a n ) > and (b\, . . . , b n ) > are the vectors of the lengths of 
the semi-axes. Then system ( 14.21) is necessary and sufficient for the inclusion S 1 ^ C Sb- 

Proof. Note first that k = I = n + 1. It is obvious that Sa Q Sb if and only if b p — a p > 
for every p — 1, . . . , n. The matrices underlying the matrix pencils A(x) and B(x) are 

A p = — (E P)n+ i + E n+ i )P ) and B p = —(E Pin+ i + E n+ \^ p ) 

dp Op 
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1 i= j = s = t, 
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(a 



67 z 



a. 





s = n + 1, j = t < n, 
i = s < n, j = t = n + 1, 
i = s <n, j = t<n, i^j, 
otherwise. 



We show that C is a solution to (14. 2p . Decompose x e R( n+1 ) 2 in blocks of length n + 1 
and write Xij for the j-th entry in the z-th block. The matrix C is positive semidefinite 
since 

n+1 



X C X ^ ^ -^i/i + 2 ^ ^ , , ^ ^ ^ ^ 3^1,1 ^n+^n+l 



i=l 



k bj 

t<j<n J 



i=l 



Oi 



^ ^ T~ + ^n+l,n+l j + 



i=l 



n 

E 

i=l 



1 - 



b 2 



xf, > 



for all x G R( n+1 ) 2 . Clearly, the sum of the diagonal blocks is the identity matrix I n+ \. 
Since every A p has only two non-zero entries, every B p is a linear combination of only two 
blocks of C, 



Bp — — Cji+i.p H Cp )n +i. 



This equality is true by the definition of C. 



□ 



Remark 4.12. Using the square matrices Eij of size (n + 1) X (n + 1) introduced in 
Section [21 the matrix C in the proof of Lemma 14.111 has the form 



£a,i <^i,2-Ei,2 

^2,1-^2,1 -^2,2 

d n ,iE n ,i 

^-C/n+1,1 ^^n+1,2 



d 



dn,n—l-E n 



n-1 



n—l,nEn—l,n 

E n ,n 
a-n TP 

h J -'n+l,n 

On 



ai jp 



a n jp 

i. 1 -'n,n+l 

On ' 



E 



n+1, n+1 



where dj 



Now we prove exactness of the criterion for the containment of a ball in a polyhedron. 

Lemma 4.13. Let Sa be a ball of radius r > in normal form (12. 4p . and let Sb be 
a polyhedron in normal form (12.31) . For the containment of Sa in Sb, system (14. 2p is 
necessary and sufficient. 

Proof. Note first that k = n + 1. Since the normal form B(x) is monic, the linear poly- 
nomials describing Sb are of the form bi(x) = 1 + J2p=i bi, P x P for i = 1, . . . , / . If Sa is 
contained in the halfspace bi(x) > 0, we have ^ > Ylp=i b\ p - 
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We give a feasible matrix C to system (14.21) to show exactness of the criterion. In this 
case, C is an (n + 1)1 X (n + 1) /-block matrix defined as follows: 

/ r 2 b 2 ■ 



1 _ rl v n 6 



JJs.t 



r b. 



rb s 







i = j < k, s = t, 
i = J ' = ft) s = t, 
i < k, j < k, i ^ j, s = t, 
i — k, j < k, s = t, 
j = k, i < k, s = t, 
otherwise. 



To show positive semidefiniteness of C, consider a vector x G IR( n+1 ^. Decompose x into 
blocks of length /, and we write x^j for the j-th entry in the z-th block. Now C is positive 
semidefinite since 



x 



s=l 



E 

i=j<k 



X 



2 a,t 



- + X 



k.s 



p=i 



+ 2 E 

i<j<k 



i<k 



E 



vi=j <fc 



\/2 



if n 



2 



p=l 



> 



for all x e R( n+1 ^. The term 1 - r 2 6 2 

is non-neeative since the ball of radius r is 
contained in Sb and therefore ^ > X)p=i & 2 P - By construction, the sum of the diagonal 
blocks is the identity matrix l[. Every B p is a linear combination of only two blocks of C, 



Bp — — C n +i )P + — C Pi n+i. 



□ 



Observe that in Lemma [4.91 14.111 and 14. 131 for rational input, C is rational as well. 



5. Block diagonalization, transitivity and containment of spectrahedra 

in polytopes 

In [16, Prop. 5.3] Helton, Klep and McCullough showed that the containment criterion 
(14. 2 p is exact in an important case, namely if Sb is the cube, given by the monic linear 
pencil 



(5.1) 



1 - 

B(x) = I 2n + - x P {Epp - E n+ p, n+p ) . 



The goal in this section is to generalize this to all polyhedra Sb given in normal form (12. 3p . 
not only for the original criterion, but also for the variations discussed in Corollaries 14.61 
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and 14.71 As in Lemma I4.9[ in case that Sb is unbounded we have to use the extended 
normal form A(x) instead of A(x). 

Theorem 5.1. Let A(x) G Sk[x] be a monic linear pencil with extended linear pencil 
A(x) G <Sfc + i[a;] as defined in equation (14.81) and let B(x) G Si[x] be the normal form of a 
polyhedron. When applied to the pencils A{x) and B(x) criterion (14. 2p is necessary and 
sufficient for the inclusion = Sa Q Sb- If Sb is a polytope then the pencil A(x) can 
be used instead of A(x). 

In order to prove this statement (where the sufficiency-parts are clear from Theorem 14 .31) 
we have to develop some auxiliary results on the behavior of the criterion with regard to 
block diagonalization and transitivity, which are also of independent interest. 

As pointed out by an anonymous referee, this theorem can also be deduced from results 
of Klep and Schweighofer in [22]. A linear scalar- valued polynomial is positive on a 
spectrahedron if and only if it is positive on the matricial version of the spectrahedron. 

We use the following statement from [TB] on the block diagonalization. As usual, for 
given matrices M 1 , . . . , M l , we denote by the direct sum j=1 M l the block matrix with 
diagonal blocks M 1 , . . . , M l and zero otherwise. 

Proposition 5.2. [TSl Proposition 4.2] Let A(x) G S k [x], B(x) G S t [x] andD q (x) G S dq [x] 
be linear pencils with D q (x) = D$ + J2p=i x pDp, q = 1, ■ ■ ■ ,m. 

If B(x) = D q (x) is the direct sum with I = Y^q=i d q , then system (14. 2 h is feasible 

if and only if for all q — 1, . . . ,m there exists a kd q x kd q -matrix C q , consisting of k x k 
blocks of size d q x d q , such that 

k 

(5-2) C q = (Cj.)f J=1 tO, Vp = 0, . . . , n : D| = £ a%C% 

is feasible. 

An analogous statement holds for criterion (14.31) and the criteria discussed in Corollar- 
ies \4-6] and \\.l\ 

Since [16] does not contain a proof of this statement, we provide a short one. 

Proof. Let C\...,C m be solutions to AO]) , and set C = 0™ 1 C q . Define C as the 
direct sum of blocks of C, Cy = ©^LxC^. Then C is a solution to (14. 2\i : C results by 
simultaneously permuting rows and columns of C and is thus positive semidefinite. We 
have B p = 0™ l Dl = 0» , ^ oyCj = £? J=1 a l3 C' y 

Conversely, let C be a solution to (14.21) . We are interested in the m diagonal submatrices 
of each block C'^ , defined as follows: For q G {1, . . . , m}, let be the d q x d q submatrix 
of C[a with row and column indices {X)r=i ^ + ^••■iEr=i^}- Now the submatrix 
C q = (Cj)fj =1 consisting of the g-th diagonal blocks of each matrix C'^ is a solution 
to (15.21) . C q is a principal submatrix of C and thus positive semidefinite. The equations 
in (15.21) are a subset of the equations in (14.21) and remain valid. □ 
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We now prove transitivity of the containment criterion. We begin with a simple auxiliary 
lemma on the Kronecker products of corresponding blocks of block matrices. 

Lemma 5.3. Let A y consist of m x m blocks of size n a x n a and B y consist of 
m x m blocks of size n& x Then Y2Tt=i(A s t ® B st ) y 0. 

Proof. First note that we have v T ^E^™ =1 A s ^ v = (v T . . . v T )A(v T . . . v T ) T > for all 

v G IR n , as in the proof of Theorem 14. 3[ and hence YlTt=i ^-st 0. 

Since A, B y 0, we have A <g> B y as well. (A st <8> B st )™ t=1 is a principal submatrix 
of this matrix and therefore also positive semidefmite. Summing up the blocks of this 
matrix and applying our initial considerations, we see that J27t=i ^st ® B st ^0. □ 

The criteria from Theorem 14.31 are transitive in the following sense. 

Theorem 5.4 (Transitivity). Let D(x) G <Sd[:r], E(x) G S e [x] and F(x) G S/[x] be linear 
pencils in n variables. If criterion (14. 2p . criterion (14. 3 p or criterion (14. 7p certifies the 
inclusion So C and i/ie inclusion Sp Q Sp, it also certifies Sp, C S^. 

The transitivity statement concerning (14. 2p can be interpreted from an operator theo- 
retic point of view. It states the well-known fact that the composition of two completely 
positive maps is again completely positive. Our approach enables us to extend the state- 
ment to the relaxed criteria (14. 3 p and (14. 7p . 

Proof. We first consider the relaxed version (14. 7p . Let C DE be the de x de- matrix certifying 
the inclusion Sp, Q Sp and C EF the ef x e/-matrix certifying the inclusion Se ^ Sp- 
C DE consists of d x d block matrices of size e x e, C EF consists of e x e block matrices of 
size / x /. 

We prove that the matrix C DF consisting of d x d blocks of size / x / and defined by 



r<DF ._ fr ,DE \ n E1 



iEF 
s,t=l 

is a solution to system (14. 3[) for the inclusion So Q Sf- 



To show C DF >z 0, we start from C DE y and C EF y 0. Define a new matrix C DE by 
{C EE )ij := (C- E ) st , permuting the rows and columns of C DE . Since rows and columns 
are permuted simultaneously, positive semidefiniteness is preserved. We think of C DE as 
having ex e blocks of size d X d. C DF now simplifies to C DF = Yll t=i ^S. E ® C EF . Using 
Lemma C DF y follows. 

Next we show F p — J2ij=i ^ifiij F b for p = 0, . . . , n. By assumption, 

(5.3) F p - J2 <fW = G EF y and E p - £ d%Cf E = G DE t 0. 

By definition of C DF and the right equation of (I5.3p . we have 

E < c r = E < E <y = e (e, - G°°) it a 

i,j=l *:i = l s,t=l s,t=l 



•iEF 
'st > 
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and then positive semidefiniteness of G EF and G DE yield 

i,j=l s,t=l 

The non-relaxed version ( 14. 2 p as well as the relaxed version ( 14. 3 p follow by choosing G EF 
and G DE in ( 15. 3 p to be zero matrices. □ 



We can now establish the proof of Theorem 15.11 which also completes the proof of 
Theorem 14.81 



Proof, (of Theorem 15.11 ) Every (monic) linear pencil B(x) in normal form (I2.3P can be 
stated as a direct sum 

i i 
B(x) = @b q (x) = 0(1/ + Bx) q . 

q=l q=l 

Therefore, Proposition 15.21 implies that system (14. 2 p is feasible if and only if the system 

k k 

C* = (C ii )} J=1 hO, 1 = E^' Vp = l,...,n: b«=J2< C ii 

i=l i , j = 1 

is feasible for all q = 1, . . . , I. Note that C q is in Sk- Hence, the system has the form 
(5.4) C q = (CV& = i hO, 1 = (h, Ci) , Vp = 1, . . . , n : W p = {A p , C«) . 

In the following, we show the existence of a solution by duality theory of semidefmite 
programming and transitivity of ( I4.2p . see Theorem 15.41 

Let b\ , . . . , 6^ be the coefficients of the linear form b g (x) = (1/ + Bx) q . Since B(x) is in 
normal form ( I2.3p . the vector b q := (b\, . . . , 6^) is an inner normal vector to the hyperplane 
b q (x) = 0. Consider the semidefmite program 

(P q ) r g := max (— b q , x) 

s.t. A(x) y 

for all q = 1, . . . , I. By assumption, (P q ) is strictly feasible and the optimal value is finite. 
Hence (see, e.g., P Thm. 2.2]), the dual problem 

min (I k ,Y q ) 

(D q ) s.t. (A p ,Y q ) = b q p Vp=l,...,n, 

Y q y 

has the same optimal value and attains it. (Note that by duality (—Ap, Y q ) = —b q .) We 
can scale the primal and dual problems simultaneously by dividing by r q and get 

l = min (I k ,Y q ) 

b q 

(D q ) s.t. (A p ,Y*) = -£ Vp=l,...,n, 

r q 

Y q y 0, 
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Yc 



in which Y q . 

Since in the dual (D q ) the optimal value is attained, in (D q ) it is as well, i.e., for all 



q 



l there exists a k x £;-matrix C q such that 
C q y 0, 1 



6 9 

(4,C«>, ^ 



As mentioned before (15.4p . the matrices C" 3 certify the inclusion C where B'(x) 
is defined as the scaled monic linear pencil 



B'(x) = 



1 



Now we have to distinguish between the two cases in the statement of the theorem. 

First consider the case where Sb is a polytope. Since B(x) is in normal form, we have 
max x£ s B (—b q ,x) = 1. Further, since Sa Q Sb, the definition of r q implies r q < 1 and 
hence, Sb 1 Q Sb- By transitivity and by exactness of the criterion for polytopes, see 
Theorem 15.41 and Lemma |4"HI respectively, there is a solution of system (14. 2p that certifies 
Sa ^ Sb- 

To prove the unbounded case in the theorem, we construct a solution to (14.21) for the 
inclusion S^ C Sjp, where B'(x) = 1 © B'(x) denotes the extended normal form (14. 8 p of 
the polyhedron Sb'- Then the claim follows by Lemma 14.91 and Theorem 15.41 as above, 
since 5*^7 C 5*^' is certified by (14.21) . 

First note that S^ C Sb is equivalent to Sa Q Sb- Denote by C the matrix that 
certifies the inclusion Sa Q Sb>- Then the symmetric (k + 1)(Z + 1) x (k + 1)(/ + l)-matrix 

k 

II II 

C := E u 












where En and the blocks 





^ 



are of size (/ + 1) x (I + 1), certifies the inclusion 



S^ ^= S-g}. Indeed, adding zero-columns and zero- rows simultaneously preserves positive 
semidefiniteness and, clearly, the sum of the diagonal blocks of C is the identity matrix 
Since in every A p the first column and the first row are the zero vector, we get 

k 



i,j=0 



0-E u + 











B 1 



v ■ 



where af, is the (i, j)-th entry of A p 



Feasibility of the relaxed criteria is again implied by the feasibility of (14.2 



□ 



6. Containment of scaled spectrahedra and inexact cases 

Contrasting the results of Sections H] and El we first consider a situation where the 
containment criterion fails and the relaxed version (14.71) is strictly stronger. In particular, 
this raises the question whether (for a spectrahedron Sa contained in a spectrahedron Sb) 
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the criterion becomes satisfied when scaling Sa by a suitable factor. In Proposition 16.21 
we answer this question in the affirmative. We then close the paper by applying this result 
on optimization versions of the containment problem. 

6.1 . Cases where the criterion fails. We review an example from [T6l Example 3.1, 
3.4] which shows that the containment criterion is not exact in general. We then contrast 
this phenomenon by showing that for this example there exists a scaling factor r for one 
of the spectrahedra so that the containment criterion is satisfied after this scaling. 

Consider the monic linear pencils A(x) — 1$ + xi(Ei )3 + E 3 i) + x 2 (E 23 + E 32 ) £ Sa[x] 
and 



B(x) 



+ Xx 



1 



■1 



+ x 2 



Clearly, both define the unit disc, that is Sa = Sb- 



Claim. The containment question Sb Q Sa is certified by criterion (14.2j) . while the 
reverse containment question Sa Q Sb is not certified by the criterion. 

First, we look into the inclusion Sb Q Sa (where the roles of A and B in (14. 2p have to 

is satisfied if and only if there exist c\ , c 2 , c 3 G K such 



be interchanged). Criterion 
that 



C 



r 1 

2 

1 
2 





1 

2 





1 

2 



1 

2 




-ci 
-c 2 



Cl 



1-C 3 



C-2 





i 

2 







1 

2 






1 

2 



D 6x6 



— Ci — c 2 
ci 1 - c 3 
L c 2 c 3 

is positive semidefinite. Since the 2 x 2-block in the top left corner is positive definite, 
the matrix C is positive semidefinite if and only if the Schur complement with respect to 
this block is positive semidefinite. One can easily check that this is the case if and only if 
ci = C3 = I and c 2 = 0. 

Conversely, Sa Q Sb is certified by (14.2D if and only if there exist ci, 
that 



, C12 6 I such 



Cl 


c 2 


c 9 


ClO 


1 

2 


c 7 








c 2 


C3 


Cll 


Cl2 


-C7 


1 

2 








Co, 


Cll 


c 4 


c 5 





c 8 








C10 


C12 


c 5 


c 6 


l-c 8 







G I 


1 

2 


-07 





1-C 8 


1 — Ci — c 4 


-c 2 - 


c 5 






c 7 


1 

2 


c 8 





-c 2 - c 5 


I-C3- 


- c 6 . 







c 



is positive semidefinite. We show the infeasibility of the system (14. 2ft . 

Assume that C is positive semidefinite. Then all principal minors are non-negative. 
Consider the principal minor 



ci 
1 
2 



1 

2 

Cl 



c 4 



Ci(l - Ci - c 4 ) 



1 

4 



Cl(l - Ck 



1 

4 



CiC 4 . 
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Since the expression in the brackets as well as the second term are always less than or 
equal to zero the minor is non-positive. Therefore, ci(l — c\) — \ = and c\c^ = 0, or 
equivalently, C\ = \ and c 4 = 0. 

Recall that whenever a diagonal element of a positive semidefinite matrix is zero, the 
corresponding row is the zero vector, that is C5 = cs = eg = en = 0. Now, we get a 
contradiction since the principal minor 

c 6 1 - c 8 

1 — Cg 1 — C\ — C4 

implies that c 6 > 2 and therefore 1 — c 3 — c 6 < — 1 — c 3 < or c 3 < — 1. This proves the 
claim. 

Now, generalizing A(x), let A r (x) be the linear pencil of the ball with radius (1 >)r > 
in normal form. With regard to the containment question Sa*- = tSa C Sb, we show the 
feasibility of system (14. 2 p for r sufficiently small. Consider the matrix 



C = 



c 








C 


r 
2 


" 





c 


— c 








r 
2 





— c 


c 








r 
2 


c 








c 


r 
2 





r 
2 








r 
2 


1 - 2c 





. 


r 
2 


r 
2 








1 - 2c . 



Obviously the equality constraints in (14. 2 1) are fulfilled. 

As above, if c = or 1 — 2c = 0, then r = 0. Therefore, < c < | and the 2 x 2-block 
in the top left corner G\\ is positive definite. Thus the matrix C is positive semidefinite 
if and only if the Schur complement with respect to G\\ is positive semidefinite. This is 
the case if and only if 

2 

1 - 2c - — > f(c) := 8c 2 - 4c + r 2 < 0. 
4c 

Assume r > \y2. Then /(c) > for all c since / has no real roots and the constant term 
/(0) = r 2 is positive. Otherwise, f(j) = —\ + r 2 < 0. Hence, system (14. 2p is feasible for 
< r < |-v/2. 

The problem of maximizing r such that the system (14. 2[) is feasible can be formulated as 
a semidefinite program. A numerical computation yields an optimal value of 0.707 ~ \y/2. 

Note that we are in the situation of Corollary 14.71 For the relaxed version ( 14. 7p . a nu- 
merical computation gives the optimal value of 0.950 ~ ||. In particular, this shows that 
the relaxed criterion ( 14. 7ft can be satisfied in cases where the non-relaxed criterion f]4. 21) 
does not certify an inclusion. 

It is an open research question to establish a quantitative relationship comparing cri- 
terion ( 14. 7p to ( 14. 2 p in the general case. 



c 6 1 
1 \ 



c 6 - 1 
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6.2. Containment of scaled spectrahedra. For a monic linear pencil A(x) G Sk[x] 
and a constant v > define 

(6.1) A v {x) := A (~) = I k + -J2x p A p , 

P =i 

the v -scaled (monic linear) pencil. Similarly, we denote by vSa := {i 6 K" : A v {x) >z 0} 
the corresponding v-scaled spectrahedron. 

Generalizing the observation from Section I6.1[ we show that for two spectrahedra Sa 
and Sb, containing the origin in their interior, there always exists some scaling factor 
v such that the criteria (14. 2p and ( 14. 3 p certify the inclusion vSa ^ Sb- This extends 
the following result of Ben-Tal and Nemirovski, who treated containment of a cube in a 
spectrahedron (in which case they can even give a bound on the scaling factor). 

Proposition 6.1. [2| Thm. 2.1] Let Sa be the cube (15. ip with edge length r > and 
consider a monic linear pencil B(x). Let fi = max p= i v .. in ranki?p. // Sa ^ Sb, then 
system ( I4.2p is feasible for the z/(/i) -scaled cube v(h)Sa, where v(n) is given by 



z/(/i) = min 




5> 



V> 2 exp (-V i du 



For all fi the bound z/(/i) > ^= holds. 

A quantitative result as presented in the last Proposition is not known for the general 
case. However, combining Proposition 16. II with our results from Sections H] and we get 
that for spectrahedra with non-empty interior, there is always a scaling factor such that 
system (14 ,2p and thus also system (14. 3p hold. 

Proposition 6.2. Let A(x) and B(x) be monic linear pencils such that Sa is bounded. 
Then there exists a constant v > such that for the scaled spectrahedron uSa the inclusion 
vSa Q Sb is certified by the systems (14. 2 j) and (14. 3p . 

We provide a proof based on the framework established in the previous sections. Alter- 
natively it can be deduced from statements about the matricial relaxation of criterion (14. 2p 
given in the work by Helton and McCullough [17] , see also [16] . Criterion ( 14. 2 p is satisfied 
for linear pencils A u (x) and B(x) if and only if the matricial version of vSa is contained 
in the matricial version of Sb- 

Proof. Denote by Sd the cube, defined by the monic linear pencil (15. ip . with the minimal 
edge length such that Sa is contained in it, which can be computed by a semidefinite 
program, see Theorem 15. II Since B(x) is monic, there is an open subset around the origin 
contained in Sb- Thus there is a scaling factor v\ > so that viSa ^ i^iSd ^ Sb- 

By Proposition 16. 1[ there exists a constant z/2 > such that for the problem V2V\Sd C 
Sb system (14. 2 p has a solution C D " B with v = v\v 2 - By Theorem 15. 1[ there is a matrix 
(j Avr > v which solves ( 14. 2 p for the problem vSa in vSd- 

Finally, Theorem 15.41 implies the feasibility of system (14. 2p with respect to vSa and 
by the matrix C A " B , as defined there. □ 
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In the proof of Proposition 16.21 we scaled the spectrahedron Sa by a certain factor 
v. Since uSa Q Sb is equivalent to Sa Q ^b, the criterion (I4.2p remains a positive 
semidefmite condition even in the presence of the factor v. Moreover, we can optimize 
for v such that the criterion remains satisfied. Proposition 16.21 implies that for bounded 
spectrahedra represented by monic linear pencils the maximization problem for v always 
has a positive optimal value. 

This yields a natural framework for the approximation of smallest enclosing spectra- 
hedra and largest enclosed spectrahedra. In [TBI Section 4], the example of computing a 
bound for the norm of the elements of a spectrahedron Sa (represented by a monic linear 
pencil) is given. This can be achieved by choosing Sb to be the ball centered at the origin, 



As we have seen in Section 16. 1[ applying criterion (14. 7p to the problem is stronger than 
specializing criterion ( 14. 2 p to it. However, for the criterion ( I4.2p . we obtain a particularly 
nice representation, it reduces to the semidefinite system 



see fl23D- 



c = (C, ^ , y o 



k 




(6.2) 



\/p= 1 



...,n, V(s,t) G {l,...,n + l} 2 : 
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